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Abstract 

The reaction 'Brd — * nAA is studied within the framework of the Faddeev 
equations as a possible tool to gain insight into the final state A-A interaction. 
The neutron differential energy spectrum gives a final state interaction that is 
sensitive to both the A-A amplitude at threshold, and the coupling between 
the A-A and 3-N channels. The latter is a result of interference between 
two mechanisms for the production of the final state, which suggests that 
this reaction could give a measure of flavor SU(3) violation in the two-baryon 
system. 
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I. INTRODUCTION 



Interest in the reaction E~d — > nAA is twofold: (i) It can be used to examine the final 
state interaction between the two A hyperons and, therefore, to gain some insight into the 
baryon-baryon interaction. The fact that this reaction requires coupling between the H- 
N and A-A channels to proceed, suggests that one may gain further understanding of the 
importance of this coupling. This reaction might possibly set a constraint on its magnitude, 
which in turn could give a measure of SU(3) violation in this system as well as the general 
two baryon system, (ii) It can be used to test the hypothesis that there is an H dibaryon, a 
six quark state. The existence of such a dibaryon would be a signature that meson-baryon 
degrees of freedom are not sufficient to describe this reaction and possibly the baryon- 
baryon system in general. If the H has a mass significantly below the A-A threshold, then 
the present reaction would provide a clean signal (a mono energetic neutron) in the neutron 
differential energy spectrum (NDES). However, if the H has a mass comparable to the A-A 
threshold, then it would be difficult to distinguish, in the NDES, between the H dibaryon 
and the final state A-A interaction. In this case one needs to compare the experimental 
data with a theoretical calculation based upon meson-baryon dynamics. To that extent this 
investigation could be used as a guide in the search for H dibaryons. 

Aerts and Dover [|I|] first examined the E~d — > nAA reaction within the context of 
estimating rates and spectra for the (H _ <f) aiom — > nH reaction. They approximated the 
decay rate for the (E~d) atom — > nAA reaction by evaluating the lowest order diagram, 
including no multiple scattering in the SNN intermediate state or between the AA pair in 
the final state. Taking into account such multiple scattering is essential if one is to model 
the enhancement in the neutron energy spectrum due to a strong AA interaction; transfer 
of the kinetic energy from the A A pair created by the E~~N — > AA transition to the neutron 
requires a complex series of multiple scatterings in the final state if the spectator neutron is 
to carry off almost all the available energy and probe the AA zero-energy scattering length. 
This is the analog of the nd — *> nnp breakup experiment (see, for example, Tornow et al. 
for a recent review) in which measurement of the proton energy spectrum was proposed as 
a means to study the nn zero-energy scattering length. The Watson-Migdal approximation 
fails to properly describe the proton energy spectrum; multiple scattering calculations ||£|] 
which include fully three-body dynamics are essential. 

To examine the reaction E~d — ► nAA within the framework of three-body dynamics, 
we must derive a set of equations for this specific problem. The new feature, unique to 
this problem, is the fact that the initial SNN system and the final AAiV system have two 
(separate) identical spin | particles, which must be in antisymmetric states. In Sec. |T| we 
first derive the expression for the breakup amplitude with both the initial two nucleons and 
the final two A hyperons in antisymmetric states. We then proceed in Sec. [II B| to adapt 
the Alt-Grassberger-Sandhas (AGS) three-body equations |5[] to generate the antisymmetric 
elastic and rearrangement amplitudes required to determine the breakup amplitude. In 
this way we minimize the number of coupled integral equations needed to extract the final 
breakup amplitude. In Sec. |I C| we write the breakup cross section as the coherent sum of 
three reaction mechanisms. To calculate the NDES for this reaction, we define the input 
two-body amplitudes in Sec. |111 A| , and demonstrate in Sec. |111 B| that the final amplitudes 
have been accurately calculated to satisfy three-body unitarity. The fact that the final state 
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interaction between the two A hyperons is part of the coupled channel AA-EN problem 
implies that there are two ways of converting the S to a A, the first is in the final state 
interaction [see diagram (a) of Fig. the second in the rearrangement amplitude [as in 
diagram (b) of Fig. |J. We find in Sec [III C| that these two amplitudes, each dominated by 
the final interaction in the AA-EN channel, are out of phase and almost cancel one another. 
It is only when the background term from diagram (c) of Fig. [j] is also included in the 
evaluation of the cross section that a relatively weak final state interaction peak is present 
in the NDES, and only if the AA-EN amplitude is dominated at low energies by a virtual 
bound state pole near the A-A threshold. The cancellation between diagrams (a) and (b) of 
Fig. [I] results in a final state interaction peak that is sensitive to the coupling between the 
A-A and the E-N channels. Finally, in Sec. [IV| we present some concluding remarks about 
the implication of the results. 



II. THEORY 



The reaction E~d — > nAA differs from the standard three-body problem, e.g. nd — > nnp, 
in that the particles in the initial state ENN differ from those in the final state A^AA, 
and in that both initial and final states contain two different pairs of identical Fermions, 
for which we need to insure antisymmetric wave functions. In addition, the EN — > AA 
conversion can take place on either nucleon in the initial state. All of these features can 
be included naturally if we work within the framework of SU(3) rather than the SU(2) of 
isospin. However, then we cannot include the mass splitting in the baryon octet, and we 
must carry out the necessary re-coupling within the framework of SU(3) algebra. To avoid 
this complication, and to take into consideration the mass splitting in the baryon octet, we 
have resorted to antisymmetrizing explicitly both initial and final states. In this section we 
first derive an expression for the antisymmetric breakup amplitude describing the reaction 
E~d — > nAA in terms of antisymmetric off-shell elastic and rearrangement amplitudes. We 
then proceed, with the help of the AGS equations to derive coupled integral equations for 
the elastic and rearrangement amplitudes. By deriving coupled integral equations for the 
antisymmetric amplitude, we reduce the number of coupled integral equations we need to 
solve. Finally, we write the cross section for this reaction in terms of these antisymmetric 
amplitudes. 



A. Antisymmetry 



As a first step in determining the antisymmetric amplitude for "Ed — > A^AA, we construct 
the initial and final states for this reaction. For the initial Ed state, we designate the two 
nucleons as particles 1 and 2 and the 5 as particle 3, while for the final state the nucleon 
can be either particle 1 or 2 depending on which nucleon the E converted into a A. We now 
introduce the antisymmetrization operator Aij defined in terms of the permutation operator 
Pij that exchanges the coordinates of particles % and j as || : 

Aij = \(l- Pij) , (1) 
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where AijAij = Aij. This allows us to define an initial state that is antisymmetric in the 
two nucleons that form the deuteron as 

\~- d) AS = A 12 \~- d) , (2) 

with the antisymmetrized state being normalized; i.e., 

AS (dZ-\~-d) AS = 1 . (3) 

We note that, if the deuteron is in the 3 Si- 3 Di channel, then An |S _ d) = |S~ d). 

For the final state we have the two A hyperons which must be in an antisymmetric state. 
However, now we need to recall that the final nucleon can be either nucleon 1 or nucleon 
2, depending on which nucleon the H converted to a A. We therefore have two possible 
configurations 0: 

|AT 1 AA) A5 = y/2A 22 \N x A 2 A 3 ) 

\AN 2 A) AS = y/ZAvilAiNtAa) . (4) 

The factor of v2 was introduced to guarantee that the antisymmetric states are normalized 
to one. Since neither of these states is the physical state, we need to take the linear combi- 
nation of these states such that the final state is antisymmetric with respect to Ni and N 2 , 
i.e. 

\NAA) AS = y/2A u \NiAA) AS 

= ^ 23 |^iA 2 A 3 ) - A 3 |Ai^ 2 A 3 ) . (5) 

In writing the above result, we have made use of the multiplication table of the permutation 
operators, and in particular the fact that 

P12 P23 — P123 — P13 P12 ■ (6) 

We will see next that the antisymmetry in the initial state between the two nucleons propa- 
gates through a symmetric breakup operator to project only the antisymmetric combination 
considered in Eq. ([5]). 

We are now in a position to define the physical amplitude for "Ed — > NAA as the matrix 
element of the breakup operator for distinguishable particles between antisymmetric and 
normalized initial and final states 0; i.e, 

T NAA ^ d = AS (NAA\U 03 \~d) AS 

= AS (NAA\ T a G U a3 \~d) AS . (7) 

a 

We note here that U 03 commutes with the permutation operator P 12 , and it is this feature 
of the breakup operator that required we write the antisymmetrized final state in the form 
given in Eq. @. In Eq. (|7|), the operators U a3 are the AGS || operators that are a solution 
of the AGS equations 

U aP (E) = SapG^iE) +J2$c n T 1 (E)G (E)U J p(E) . (8) 
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Here 5 a p = 1 — 8 a p, while Go is the free three-body Green's function, and T a the amplitude 
for particles ft and 7. 

We are now in a position to write the physical amplitude for 3d — > NAA in terms of a 
linear combination of the AGS operators U a p to maintain the antisymmetry in both initial 
and final states. We then can use the AGS equations to derive integral equations for these 
antisymmetric combinations of U a p. The advantage of this procedure is a reduction in the 
number of coupled integral equations we need to solve to construct the breakup amplitude. 

Combining Eq. (|5]) with Eq. ([?]), we can write the physical amplitude as: 

TjVAA^Hd = (N 1 A 2 A 3 \A23T, T ^ G oU a 3A 12 \Sd) 

a 

- (A 1 N 2 A 3 \A 13 J2T a G U a3 A l2 \Ed) . (9) 

a 

The permutation operator P\ 2 exchanges the coordinate of particles 1 and 2. In the final 
state, this interchanges the position of the particle in the ket and the label on the nucleon; 
i.e., 

(A 1 N 2 A 3 \ = (A^iA 2 A 3 | P12 ■ (10) 

With the help of the identities 

P12 A 13 = A 23 P12 PnA 12 = -A 12 (11) 

that follow from the multiplication table of the permutation operators, we can write the 
breakup amplitude as 

T NA A<-nd = 2 (JViA 2 A 3 | A 23 A 12 U 03 A 12 \Ed) 

= 2(N 1 A 2 A 3 \A 23 A 12 T a G U a3 A 12 \Zd) . (12) 

a 

In writing the above expression for the breakup amplitude, we have maintained the antisym- 
metry operators on both sides of the breakup amplitude U 03 . To reduce the right hand side 
of the above expression for the breakup amplitude in terms of antisymmetric two-cluster 
final state amplitudes, we make use of the fact that 

T a P a p = P af3 T/3 , and T 7 P af3 = P afi T 7 for a ^ 7 ^ ft , (13) 

and operate with A\ 2 in the final state on the breakup amplitude. This will assure us that 
the final state is antisymmetric in particles 1 and 2. The resultant breakup amplitude is the 
sum of two terms. The first has a final state interaction (FSI) between the two A hyperons, 
while the second term has an A^A FSI; i.e., 

T NAA ^ d = (A^A 2 A 3 | A 23 T x G [ U 13 - P 12 U 23 } A 12 \Ed) 

+(^ 1 A 2 A 3 |^ 23 T 2 Go [ (U 23 - P 23 U 33 ) - P 12 (U 13 - P 13 U 33 ) ] A l2 \Zd) . (14) 

The linear combination of the amplitudes in Eq. ([14]) insures that the antisymmetry in the 
final state is preserved at the operator level. 
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The AGS equations take the form of a closed set of coupled integral equations for physical 
amplitudes only if the two-body interaction is assumed to be separable. Since we will be 
using separable interactions for calculating the cross section for this reaction, it is most 
convenient to introduce this approximation at this time as it will give a physical meaning 



to the terms resulting from Eq. Production of the AAA final state will require the 

conversion of HA into AA. This can take place either in the two-body amplitude T\ in the 
FSI where the conversion takes place on nucleon 2, or in the three-body AGS amplitudes 
[U\3 — P\2 ^23]- To expose the mechanism for conversion, we need to write the amplitude 
Ti in the form of 2 x 2 matrix. In three-body Hilbert space this takes the form 

T 1 (E) = Y: \g c a1 N l )r c aP {E-e 1 )(gl- 1 N 1 \, (15) 

where the sum a, (3 runs over the coupled channels AS, AA, and the superscript c indicates 
that this is a coupled channel partial wave. Since we have written the two-body amplitude 
Ti(E) in three-body Hilbert space, e\ is the energy of the spectator nucleon 1. 

Because the AA is an isospin zero system, this matrix structure for 7\ (E) is only present 
for this isospin channel and partial waves in which the space-spin wave function are anti- 
symmetric. In all other partial waves the HA system does not couple to the AA channel, 
and the corresponding amplitude is a single channel amplitude, and has the same form as 
the AA interaction, which can be written as 

Ti{E) = [gNA^Ni) t na {E - ti) (g NAj : JV,| j + 1 = 2, 3. (16) 

In writing the separable representation for the AA amplitude, we have excluded the particle 
label from the tn A (E) because this quantity is the same for the nucleon interacting with 
either A. 

Making use of the above separable representation for the two-body amplitude, we can 
write the amplitude for AAA <— "Ed as 



Tnak^n = (N 1 A 2 A 3 \ A 23 \g c AA ; 



T AA;SjV X N:S + T AA;AA Yn-E 



+ 2 (N 1 A 2 A 3 \ A 23 \g N A] A 2 > r NA Y A . 3 . (17) 

The first two terms on the right hand side correspond to the final interaction occuring 
in the AS-AA channel, while the last term has the interaction in the AA channel. The 
antisymmetric AGS amplitudes X a p and Y a p are defined as 

Aat ; = = — j= (g~N] Ai| Go [Ui 3 — P12 U23] G \gNN', ^) 

= 7/2 [ Xni ' s ~ XN r,s) (18) 

for the antisymmetric (HA) A <— Ed, while the antisymmetric amplitude for (A A) A <— Ed 
is defined as 

Yn-,~ = —7= (#aa; Ai| G [U 13 - P12 U23] G \g N N] S) 



V2 
1 
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[Y Nr , s -Y N2 . 3 ] . (19) 
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In writing the above two definitions for the amplitudes, we have dropped the superscript as 
the definitions are valid for both coupled and uncoupled channels. Finally the antisymmetric 
amplitude for the reaction (A/A) A <— Ed is given by 



>A;s = - (^A 3 ; A a | Go [ (U 23 - P 23 C/33) - P12 (U 13 - P 13 Uaa)] G \g NN ; E) 



1 
2 



y A 2 ;S 



r A 3 ;S 



V N 2 1 



(20) 



In defining X a p and Y a p, we have made use of the fact that in a separable approximation 
we can write the initial state \Ed) = Go \gd] E). In addition, it is assumed that the deuteron 
is in an antisymmetric state; i.e., Au\gd', S) = \gd] E). The amplitude given in Eq. (pp can 
be represented diagrammatically as in Fig. |l[ 



B. The AGS equations for Ed -> AAA 

We are now in a position to derive integral equations for the antisymmetric amplitudes 
X a/ 3 and Y a p defined in Eqs. (|i"8D -(pO"D, and required in Eq. (|i~7|) to construct the total 



amplitude for breakup. In addition to these amplitudes, we need the amplitude for Ed 
elastic scattering, which is basically a matrix element of U 33 . Thus for X^.^ we have, after 
making use of the AGS equation, 

Xpj-E = ~~fK (9sn'i Ai| Go [U13 — P\2 U23 ] Gq \gNN', E) 



V2 
1 

1 
1 
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— (9l 3 iv 2 ;Ai|(l-Pi2)G | 5ArA r;H) 

+4; telstfa ; Ail Go T 2 Go ( C/23 - Pia ^13 ) G 1^; S) 
+4= (g c S3N2 ;N 1 \(l-P 12 )GoT 3 G U 33 Go\g NN ;E) . (21) 



Introducing the separable representation for the amplitudes T 2 and T 3 allows us to turn the 
above expression into an integral equation of the form 

Ajv ; s = ^jV;H + Z N;N T HW Atv ; S + ^N;N T SN;SN -^N;S 

+ Z N C ,N T §.N;AA Yn;S + Z N;3 T NN A H;S , (22) 

where the antisymmetric elastic amplitude is given by 

As ; s = (gNN', ^1 G ^33 G \gNN', ^) , (23) 
while the antisymmetric Born amplitudes Z a p needed in Eq. (j2~2"D are defined as 

Zn-,E = — ^ (9Z 3 N 2 '■ Ai I (1 — Pi 2 ) G Ifi'ATAr; H) 

= ^ [^;H-^iV 2 ;s] (24) 



and 

Zn-,n = —(gs s N 2 i Ni\ G Q \gsN] N 2 ) 

= _ Z Nl . N2 . (25) 

In Eq. (p2j) we have the amplitude X^. s in terms of the amplitudes X^. s , Y$. s , Xn ]S , 
and To close this set of integral equations, we generate an equation for each of these 

amplitudes following the same procedure we adopted for deriving Eq. (p2|). For the amplitude 
Xn-,z, we have by definition the same expression as in Eq. (|2ll) without the superscript c . 
This results in an equation for the Xjf-s amplitude which is identical to Eq. (p2"|) without 
the left hand superscript c . 

For the amplitude Y^. E , we have 

^JV;H = 7^ (fl'AA) ^"ll Co [ ^13 — -^12 ^23] G \gNN', E) 

(fl'AA) ^ll ^0 ^2 G*0 [ ^23 — -Pl2 Ui 3 ] G \gNN', 2) 



v/2 
V2 



+4^ (gl^N^GoTzGo [U 33 -P 12 U 33 ] G \g NN ;E) . (26) 



Making use of the fact that T 3 = P 23 T 2 P 23 and {g^ A ; Ni\P 23 = —(g^N^, the integral 
equation for Y N]S reduces to 

Yn-,~ = Zn-a t na *a ; e , (27) 

where 

Z N . A = V2 (g c A \; Nx\ G \g NA ; A) 

= v/2 Z% vM . (28) 

In this same way we establish an equation for Y A . S to be 

1 A;S — ^A;N T AA;~N ^N;S ~>~ ^A;^ r AA;AA 1 N;S 

+Z A . A r NA Y A . E , (29) 



where Z A . N is defined as 



Za-,n = (0iViAa; A 3 | (1 — P 23 ) G |#aa; Nx) 

= y/2 Z A2]Nl , (30) 



and 



J A;A — ~ Z/ A 3 ;A 2 



(31) 



Finally to close the set of coupled equations, we must write the equation for the elastic 
amplitude. Making use of the definition of Xs-s, given in Eq. (|23|) , we get 



S 



where 



Xe-,E — Z-z-n T-en Xn-E + Z~^ N 7"|jV;HAf -^iV;S 

+ ^s';AT T EN;AA ^N;E J (32) 



Ze-,n = ^= (9nn, S| (1 - P 12 ) G \gEN] N) 

= [ Zez-Nx - Ze 3 ;n 2 ] • (33) 



Equations (^2|), (|27|), (p9[) , and ( |33| ) now form a closed set for the two-cluster to two-cluster 
antisymmetrized amplitudes required to construct the breakup amplitude and, therefore, 
the cross section. In the appendix we provide a summary of the antisymmetrized Born 
amplitudes in terms of the single particle exchanged amplitudes. 



C. The cross section for Ed — * NAA 

We now turn to the determination of the cross section, and in particular the neutron 
energy spectrum for the capture of S~ on the deuteron in terms of the amplitudes defined in 
Sec. [11 A|. Since our ultimate aim is to gain some insight into the low energy AA interaction 



in a reaction in which the S~ is captured from an atomic orbit, we will restrict our analysis 
to S-wave two-body interactions. With this simplification, the only angular momentum in 
the problem is the spin of the three particles, and the total spin will be either | (doublet) 
or | (quartet). However, because the two A hyperons are in a 1 S state, only the doublet 
channel contributes. In addition, by including isospin and taking into consideration the fact 
that the two A hyperons and the deuteron are in an isospin zero state, the total isospin is 
| with a projection — |, corresponding to a S~ in the initial state. With these restrictions, 
we can write the breakup amplitude as the sum of the three diagrams in Fig. [I] in the spin 
isospin | state as: 

T A An^E-d = Yl {\m S2 \m Si \S N Ms N ){S N Ms N \m Sl \\M s ) 

S N M s =±i 

x [ A eN .. iB + B £n . /s + C lN , ts ] (S d M Sd \m as \\ M s ) , (34) 

where £ a = {S a ,t a }, with S a and t a the spin and isospin of the pair (^7) respectively. 
This basically defines the channels in the two-body subsystems. The amplitude Ae N -e s 
corresponding to Fig. |l|(a) is given in terms of the AGS amplitude X a/3 by 

A^ ;fe ( Pl ,k i; k°) = g iN { Px ) t aa .,en[E - 6(h)] X^k^kl) , (35) 

V Z7T 



where k| and ki are the initial on-shell momentum of the S~ and the final momentum of 
the neutron respectively. The momentum pi in Eq. (|35| ) is the relative momentum of the 
two A hyperons in the final state. The AGS amplitude Xi N] i s is that defined in Eq. ( |T%D for 
a state with total spin and isospin |. Here we note that the £5^0 restricts the contribution 
of Fig. |I](a) to the amplitude in which the AA interaction occurs in the 1 Sq state. The 
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amplitude corresponding to diagram Fig. |l|(b), in which the conversion from SiV to AA 
takes place in the multiple scattering before the FSI, is given by 

B iN] UPi^i-y s ) = %= 9iM taaME - e(h)} ^(k^kl) , (36) 
V In 

where the AGS amplitude Yt Ni i s is defined in Eq. (P^|) and is non-zero in total spin, isospin 
| with the two A hyperons in the 1 So state. If we compare the full breakup amplitude given 
in Eq. ( |TTD and (0), we observe that the antisymmetrization operator A23 in the first two 
terms on the right hand side of Eq. ([17]) are not present in the expressions for A^ N .^ and 
Bt N .g s . This is due to the fact that the AA interaction in the final state has been restricted to 
those channels that satisfy the Pauli principle for the A hyperons. This was simple to achieve 
by proper choice of partial waves. Unfortunately, that is not possible for the last term in 
Eq. fllTD , and in this case we need to include the antisymmetrization operator explicitly. In 
addition, we need to maintain the same spin coupling as the diagrams in Figs. 0(a) and 0(b). 
This will involve a re-coupling of the spin, with the resultant amplitude given by: 



c^ s ( Pl ,k i; ]4) = £ (-if s N & J s * DN 



'A II 
5 A I 2 2 



1 1 c 

Sa 



X 



[Ci2 fe (p 2 ,k 2 ;k°) + (-1)* Cg] fe ( P 3,k 3 ;k 3 )] (37) 
where a = \/1a + 1, R — Sa and, R' — 1 — Sn- The phase (— 1) R ' ensures that the two 



A hyperons are in states that satisfy the Pauli principle. The amplitude Cl*.« , i = 2,3 is 



given by 



C« fc (p,, k,; k°) = -±= g lA { Pi ) r AN [E - e{h)} Y th ,Uh kl) , (38) 

- '47r 



where the AGS amplitude Y^ A .£ B is defined in Eq. (|20|). 

To calculate the neutron differential energy spectrum (NDES) we square the breakup 
amplitude, perform spin averaging over the incoming particles, and integrate over all the 
kinematical allowed final states for the two A hyperons. This gives the cross section as a 
function of the final neutron energy. In the three-body center of mass system, the NDES 
takes the form 

d 3 a (27r) 4 /Z2d&i7Bjv 



d 3 k 2 d 3 k 3 6(kx + k 2 + k 3 

2 U2 



x ^-^-^-^ )|M(kl ' k2 - k3;k|)|2 - (39) 

where the (^-functions maintain energy momentum conservation. The amplitude \M\ 2 is 
given as the square of the sum of the diagrams in Fig. |1|; i.e., one obtains 

\M\ 2 = W \A tN]ls + B tN . ts + G A; , S | 2 , (40) 

6 Si 

with Ae N -e s , Be N -e s , and Ci A] e s given in Eqs. (|35D , d36|) , and d37|) . Because we are summing 
over the spin projections and integrating over the momenta of the two identical A hyperons, 
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we have included a factor of 1/2! in Eq. (0). This factor is physically crucial as it avoids 
double counting of states, that afterall represent only interchanges of identical particles. The 
integral in Eq. ([JO) can be reduced to a two-dimensional one which needs to be evaluated 
with some degree of accuracy &[L1 



III. NUMERICAL RESULTS 

Having established that the neutron differential energy spectrum (NDES) is basically 
the square of the coherent sum of the diagrams in Fig. [I], we are now in a position to 
determine whether the reaction S~d — > nAA is, in fact, a means to investigate the final 
state AA interaction. To achieve this we must: (i) Solve the coupled integral equations 
for the amplitudes X N , S , Y N;3 , and Y A , E , i.e. Eqs. (0), (|27|), (|2|), and (|32]), and then 
construct the breakup amplitudes Ai N g 3 , Bg N .£ s , and Ci N] i s defined in Eqs. (^), ([36D, and 
(|T7|). (ii) Establish the existence of an enhancement in the FSI region that is sensitive to 
the choice of the AA interaction. 



A. The input amplitudes 

The proposed reaction is initiated by H~ capture from an atomic s orbit; as a result 
we consider a low energy reaction in which the two-body interactions are dominated by the 
S'-wave contribution. Inclusion of the Coulomb interaction in the initial state to get the S~ 
in an atomic orbital turns a relatively simple three-body problem into a much more difficult 
numerical calculation due to the complex analytic structure of the Coulomb amplitude. 
Because we are predominantly interested in the final A-A interaction, we have assumed 
that the initial S~ is incident on the deuteron at an energy of one MeV, which is below 
the deuteron breakup threshold. This allows us to neglect the Coulomb potential, since we 
don't have an initial bound state, and will enable us to check the numerical accuracy of our 
amplitudes with the help of three-body unitarity. 

Since we are assuming an 5-wave interaction and considering an initial state below the 
breakup threshold for SNN, we may assume that all orbital angular momenta in the problem 
can be set to zero. As a result, the partial wave expansion for the elastic and rearrangement 
amplitudes reduces to 

X la . e0 (k a ; kg) = — X £a . i0 (k a ; hp) . (41) 

With this limitation on the angular momentum, the coupled integral equations can be re- 
duced to a set of one dimensional equations with the channel £ a labeled by the total spin 
and isospin (S a ,t a ) of the pair 7). In Table | we present the nine channels included 
in the solution of our coupled integral equations for total spins and isospin |. Because 
the kernel of these equations have the standard moving singularities encountered in any 
three-body breakup problem, we have followed the same procedure implemented in the n-d 
breakup reaction using the rotation-of-contour method [|13|,[|,|ll[] . To test the accuracy of 



our numerical procedures in solving the coupled integral equations and the construction of 
the breakup amplitudes corresponding to the three diagrams in Fig. [I], we have made use of 
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three-body unitarity to compare the imaginary part of the elastic E-d amplitude to the total 
cross section. To calculate the total cross section for elastic Ed scattering, and in particular 
the imaginary part of the elastic amplitude, we also need to solve the integral equation for 
the spin quartet, or S = | case. In Table || we have the three-body channels for the S = | 
amplitudes. 

To evaluate the NDES, we must define the input two-body separable potential. From 
Tables | and [TJ it is clear that we need the two-body interactions in the following channels: 



NN 
EN 
IN — AA 



1 C 3n 
Do, 0\ 

1 S t = C 



i S 1 t = 0, 1 



(42) 



For the NN potentials we have used the S'-wave separable potentials of Afnan and Gibson || , 
while the EN and EA^-AA potentials are those constructed to give the same effective range 
parameters as the SU(3) rotated Nijmegen model D [1C] with soft core fli"l"| , |i"2"] . Since we 
will examine the sensitivity of the FSI peak to the choice of the AA-SA" interaction, we 
present in Table [TIT] the effective range parameters and the A-A binding energy for the four 
potentials under consideration. Here we note that potentials SA and SB do not support a 
A-A bound state, while potential SCI and SC2 do support a bound state. In particular, 
the potential SB has a A-A scattering length that is comparable to the 1 So n-n scattering 
length; i.e., it has a virtual state with a pole in the amplitude on the second energy sheet 
near the A-A elastic threshold. 



B. Three-body unitarity 

Three-body unitarity for this reaction gives a relation between the imaginary part of the 
elastic Ed amplitude and the total cross section at the corresponding energy; i. e, one has 

Im[X E . s ] = -— — v Ie \ OTot , (43) 

where the relative velocity t> re i = k%/ ^sd, and the total cross section is the sum of the elastic, 
breakup, and reaction cross sections, or 

0"Tot = Oel + OB-up + CW . (44) 

We can now determine the total cross section o"Tot from unitarity, while the elastic scattering 
cross section can be determined by integrating the elastic differential cross section. This gives 
us, in the absence of a A-A bound state, the total breakup cross section as 

0"B-up = 0"Tot - Cel , (45) 

while in the presence of a A-A bound state 

°"B-up = OTot — CT e i — CT rcac , (46) 

with cr rcac being calculated by integrating the rearrangement differential cross section for 
E + d -> n + (AA). 
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As a test of the numerical accuracy in the solution of the integral equation, we compare 
the results for the total breakup cross section as calculated from unitarity and the elastic 
and reaction cross section (i.e. Eq. (f45|) or Eq. (^) ), with the result of integrating the 
NDES over the neutron energy. In Table |V| we compare the results for these two methods 
of calculating the total breakup cross section. We have included results for two of the four 
potentials we will consider. Potential SB has no A-A bound state, while potential SCI has 
a bound state, and therefore a rearrangement cross section. The comparison is done at an 
energy of 25.7 MeV above the iVAA threshold which is just below the SNN threshold. The 
good agreement between the two methods used to construct the total breakup cross section 
is a clear indication that the numerical procedure for calculating NDES is satisfactory. This 
is particularly ture when the unitarity approach for calculating o"b-u P involves taking the 
difference between the total and total elastic cross sections with both being an order of 
magnitude larger than the total breakup cross section. 



C. The neutron differential energy spectrum 

We are now in a position to examine the NDES for the different interactions. In Fig. |2] we 
present the neutron differential energy spectrum for the four potentials under consideration. 
Here we observe that for low energy neutrons all four potentials give basically the same 
shape for the cross section with the broad peak for low energy neutrons being independent 
of the choice of potential for the AASN 1 So channel. On the other hand in the FSI region 
(i.e. for large neutron energies) the four potentials give different results with potential SB 
having a small FSI peak. Considering the similarity between the potential SB and the 1 So 
n-n potential, it is surprising that we don't see a large FSI peak similar to that observed in 
the n-d breakup reaction. In the n-d breakup, the enhancement in the cross section at the 
end of the proton spectrum is a result of the fact that the energy available for the final n-n 
interaction is small, and therefore the cross section is dominated by the pole in the two-body 
n-n amplitude near zero energy which is known as the 1 So virtual state. The questions then 
are: (i) Why don't we see a similar FSI peak in the present reaction? (ii) Why is the broad 
peak at the low energy part of the spectrum the same for the four potentials? 

To answer the first question, we consider the cross section for potential SB in more 
detail. Since the FSI peak should be the result of a pole near zero energy in the AA-5A" 
1 S amplitude, we expect a contribution to this FSI to come from diagrams (a) and (b) of 
Fig. [TJ We therefore consider the NDES for the following combination of diagrams: (i) The 
diagrams (a) and (c) of Fig. [I], see Fig. [| (ii) The diagrams (b) and (c) of Fig. [I], see Fig. £|. 
Here we observe that in both cases there is an enhancement in the FSI as we originally 
expected. However the magnitude of the FSI peak is considerably larger in Fig. |] than in 
Fig. |3|. This suggests that we have constructive interference between diagrams (b) and (c), 
and destructive interference between diagrams (a) and (c) of Fig. [I]. In this way we establish 
the fact that diagrams (a) and (b) of Fig. [I] are of opposite sign leading to a cancellation 
between the two breakup amplitudes that provide the enhancement due to the final state 
A-A interaction. 

To establish the relative contribution of diagrams (a) and (b) of Fig. [TJ, we present 
in Figs. [5] and |6| the NDES for diagrams (a) and (b) respectively. A comparison of the 
magnitude of the FSI peak in these two diagrams establishes the fact that the contribution 
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of these two diagrams to the total amplitude is the same. This is expected considering the 
fact that it is the same pole that gives the FSI peak in the two figures. From the above 
analysis we may conclude that diagrams (a) and (b) almost cancel each other and the only 
reason we have any FSI peak in the total amplitude is a result of the fact that diagram (c) 
enhances diagram (b) relative to diagram (a). The obvious difference between this reaction 
and that for the n-d breakup is the fact that in this reaction we have a coupled channel 
problem for the FSI (i.e., AA-EN), while in the n-d case we have only one channel for the 
final n-n interaction. Finally, we note that this cancellation does not depend on the sign 
of the coupling Ca ; h between the AA and EN in the potential, but requires that the same 
interaction be used for the FSI as is used in generating the inelastic and breakup amplitudes. 
To understand this we note that a change in the sign of Ca ; h changes the sign of t AA; zn 
but not the sign of taa ; aa, and since each term in the multiple scattering series for each of 
the three amplitude has an odd number of t aa . sn , the sign of all three amplitudes changes 
under the substitution 

CA;H — ► — CA;H • 

As a result the cross section does not give us any insight into the sign of the coupling between 
the AA and EN channels. 

We now turn to the broad peak in the NDES at the low end of the neutron spectrum. 
Here, we first observe from Fig. |6] that diagram (b) does not contribute to this peak. The 
fact that it is not present in the NDES resulting from diagrams (b) and (c), see Fig. |4], 
suggests that the main contribution to the broad peak at the low energy end of the neutron 
spectrum is due to diagram (a), which in lowest order is given by 

9n(Pi)t aa . sn X n ^ « gN(Pi)T AA . SN Z NjE , 

and corresponds to the E~ interacting with the proton in the deuteron to generate two A 
hyperons. In this case the neutron spectrum is determined by the momentum distribution 
in the deuteron. To illustrate this we present in Fig. [7| the lowest order contribution from 
diagram (a) of Fig. p] to the NDES in which the FSI is almost non-existent. However, as 
we add the higher order contributions in the (divergent) multiple scattering for Xjv ; h to the 
NDES, the low energy peak is reduced in magnitude and is converted from a sharp peak to 
the broad peak we found in the full calculation. From this we may conclude that the broad 
peak in the neutron spectrum is to a large extent determined by the momentum distribution 
in the deuteron, and since the same deuteron wave function was used to generate the results 
for the four AA-SA" potentials in Fig. 0, the shape of the peak is basically the same. The 
difference in magnitude is due to the difference in the AA-HA" potentials used. 

The fact that the FSI peak is suppressed as a result of the cancellation between the 
diagrams (a) and (b) of Fig. [I], suggests that the height of the peak might be sensitive to the 
coupling between the A-A and E-N channels, since the absence of such coupling in the FSI 
would have given us the result in Fig. ^ in which the NDES is dominated by the FSI peak. 
The potential SB is based on a separable approximation to a one boson exchange potential 
(OBE) in which the long range part was determined by the SU(3) rotation of Nijmegen 
model D potential. The short range part of the interaction was adjusted in such a way that 
the SU(3) rotated OBE potential was not altered outside of 0.8 fm ||IT| , |T2|| . Because the 
long range part of the coupling potential is determined by the K- and A'*-exchanges which 
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in turn is determined by £77(3) except for the masses of the mesons and baryons, we expect 
the degree of cancellation between diagrams (a) and (b) may be determined by £77(3). In 
an attempt to test this idea, we scaled the coupling between the A-A and E-N channels; 
i.e., we assumed 

Cah — ► R Cas , 

where Cah is the strength of the coupling between the A-A, and the H-7V channels in the 
separable potential. Here, we chose the scale factor R = 0.5, 0.75, 1.0, and 1.25. In this 
way we could investigate changes in the FSI peak with variation in the coupling between 
the A-A channel and the S-iV channel. To maintain the position of the virtual pole in the 
amplitudes fixed, we adjusted the other parameters of the potential to retain a scattering 
length oaa ~ 21 fm. 

In Fig. ^ we present the NDES in the FSI region for different values of R. We note 
here that the R = 1 case gives the largest FSI peak, with the other values of R invariably 
reducing the magnitude of the FSI peak. This suggests that a measurement of the NDES 
could give more than the A-A scattering length in that the height of the peak could shed 
some light on the strength of the coupling between the A-A and S-7V channels and, in this 
way, could test the degree of violation of £77(3) symmetry in the two baryon system. If the 



suggestion of Dover and Feshbach |Tj] that there is much less dynamic SU(3) symmetry 
breaking in the AA — > SiV than in AA — > AA, then the height of the FSI peak might be 
used to explore this question. 



IV. CONCLUSIONS 

The main motivation for examining the reaction — ► nAA was to determine whether 
this reaction might provide some insight into the A-A interaction at low energies and possibly 
give a measure of the scattering length. In this way we can establish whether the A-A and 
the n-n belong to the same flavor multiplet as predicted by £77(3). At the same time if 
this reaction is to be used to establish the existence of the H dibaryon, then we should 
establish some quantitative measure of what the cross section would be in the absence of the 
H. In particular, we need to establish which features of the cross section can be reproduced 
without the introduction of the H particle. 

From the results of the above calculation we have established that the neutron differential 
energy spectrum does, in fact, give a FSI peak in the event that the A-A amplitude has 
a pole on the second Riemann sheet of the energy plane near the A-A threshold. Thus, 
experimental observation of a FSI peak is a strong signature that the A-A interaction is, 
in fact, similar to the n-n interaction as suggested by flavor SU(3). In contrast to the 
n-n interaction, the A-A interaction is part of the coupled channel AA-SiV, and we have 
established that in the case when the FSI is a coupled channel problem, then there are two 
distinct diagrams that contribute to the FSI peak. In this case the two diagrams (a) and 
(b) of Fig. H almost cancel each other resulting in a suppression of the FSI peak. In fact the 
background contribution from the diagram (c) in Fig. [TJ effectively preserves the final state 
interaction by interfering constructively with diagram (b) and destructively with digram (a). 
This cancellation is a unique feature of the coupled channel nature of the FSI. 
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The presence of the interference between diagrams (a) and (b) of Fig. [I] makes this 
reaction a possible tool to measure the strength of the coupling between the A-A and E-N 
channels. Since this coupling is determined in the OBE model by the K- and A^-exchanges 
whose coupling to the baryon is determined by SU(3), we expect the magnitude of the FSI 
peak to be a candidate for the determination of SU(3) violation in the two baryon system. 

Finally, we observe that the low energy part of the NDES has a broad peak which is 
largely determined by the momentum distribution in the target deuteron. Although this is 
not the ideal reaction to examine the deuteron wave function, the fact that we have a major 
cancellation between the three diagrams that contribute to the FSI suggests that the FSI 
peak might be sensitive to the choice of deuteron wave function. In the present investigation 
we have used the simplest possible 3 Si deuteron, and a more realistic wave function for the 
deuteron (e.g. one with short range repulsion) could give a different result for the magnitude 
of the FSI peak. The magnitude of the FSI peak could also be sensitive to the choice of 
the E-N and the AASN interactions. For the present investigation we have chosen to 
use simple separable potentials that give the same effective range parameters as the OBE 
potentials. The existence of experimental data would be required to justify the extension of 
the present results to include more realistic two-body input. 
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APPENDIX A: THE KERNEL OF THE AGS EQUATIONS 

To solve the coupled integral equations for the elastic and rearrangement amplitudes 
Xg a .i g and Yi a .£ , we need to define the kernel or the Born amplitudes Zi a .^ p . Since the 
reaction under consideration involves two Hilbert spaces, i.e the SNN and AAN, and the 
Born term does not couple these two spaces, we can define the Born term for each space 
separately. Thus, in the SNN Hilbert space we have the nucleon exchange amplitude Z£ N] e s 
and the H exchange amplitude Ze m e N - To evaluate these amplitudes we need to relate them 
to the corresponding amplitudes in which the the initial and final state are related by the 
cyclic permutation of the particle label for which we have a standard expression ||. Thus 
the A^-exchange is given by: 

Ze N -e s = [Zi a -i N ]^ = —j= [Z Nl . 3a - Z N3 . Sa ] 

= V2Z 1;3 = V2((23)l|Go|(12)3> , (Al) 

where particles 1 and 2 are the nucleon, while particle 3 is the S. The S exchange Born 
amplitude can be written in a similar manner as: 
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Zl N ;£ N = —Z Nl - N2 — ( — 1) + Zl-2 

= (-l) p+1 ((23)1|G |(31)2> . (A2) 

Here the phase P results from the exchange of the coordinates of particles 1 and 3 in the 
coupling coefficient and is given by P = si + s 3 — S 2 + T\ + r 3 — t 2 , where Sj and n are the 
spin and isospin of particle i, which in this case is ~, while S a and t a are the total spin and 
isospin of the pair (Pj). 

On the other hand, in the AAiV Hilbert space, we need the A^-exchange Born amplitude 
Zg N -e A and the A exchange amplitude Zi A -£ A . These are given in terms of the cyclicly defined 
Born amplitudes as: 

Ze N] e A = v2ZjVi;A 2 = V%Zi ; 2 

= y/2 ((23)1|G |(31)2) , (A3) 

and 

Z £a;£a = -^a 3 ;A 2 = (-!) P+1 Z 2-,3 

= ((31)2|G |(12)3) . (A4) 

In this case the the phase P' results from the exchange of the labels on particles 1 and 2, 
and is given by P' = si + s 2 — S 3 + T\ + r 2 — £3. The cyclic order Born amplitudes can in 
general be partial wave expanded for any three particles with specific spin and isospin HQ . 
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TABLES 



TABLE I. Three-body channels allowed for the S = ^, T = configuration. 







Channel 




t a 


da 


E(NN): 




NN ( Si) 
NN (%) 





1 


i 




N(EN): 




EN (%) 
EN ( 3 Si) 
EN ( 3 Si) 




I— 1 O K o 


o 



1 
1 


JV(AA): 




AA OSo) 










A(NA): 




AN (%) 
AJV ( 3 Si) 




l 

2 
1 
2 




1 




TABLE II. Three-body channels allowed for the S = 


|, T = i configuration. 






Channel 






S a 


E(NN): 




NN ( 3 Si) 







1 


N(EN): 




EN ( 3 Si) 
EN ( 3 Si) 





1 


1 
1 


TABLE III. The effective 


range parameters for the AA-EN 


coupled channel 


separable poten- 


tials in the 1 S'o partial wave. 










Pot. 


<2aa (fm) 


taa (fm) 


azN (fm) 


r~ N (fm) 


B.E. (MeV) 


SA 


-1.90 


3.33 


-2.08-0.81i 


3.44-0.22i 


UB 


SB 


-21.0 


2.54 


-2.07-6.52i 


2.62-0. 15i 


UB 


SCI 


7.84 


1.48 


3.05-5.28i 


1.45+0.074i 


0.71 


SC2 


3.36 


1.0 


3.35-2.50i 


1.83-0.10i 


4.73 
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TABLE IV. The total breakup cross section as calculated via unitarity (o"b_ up ) and by inte- 
grating the NDES (a^_ ) at an energy of 25.7 MeV. 



AA - EN 


°B-up 


a B-up 


SB 


83.32 


83.77 


SCI 


111.69 


111.69 
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FIGURES 





(a) 



(b) 



(c) 



FIG. 1. The three distinct contributions to the breakup amplitude, with (a) and (b) corre- 
sponding to the first term in Eq. ( |l7| ) and (c) corresponding to the last two terms in Eq. (|l7|). 




FIG. 2. The NDES following the capture of 3 on the deuteron for the potentials SA, SB, SCI, 
and SC2. 
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FIG. 4. The NDES following the capture of S on the deuteron for the potential SB including 
diagrams (b) and (c) of Fig. [[]. 
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FIG. 5. The NDES following the capture of on the deuteron for the potentials SB including 
diagram (a) of Fig. |l]. 
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FIG. 6. The NDES following the capture of H on the deuteron for the potentials SB including 
diagram (b) of Fig. |l[ 
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FIG. 7. The NDES following the capture of E on the deuteron for the potentials SB including 
diagram (a) of Fig. |l] to lowest order, i.e TAA ; sjv ^iV;S- 
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FIG. 8. The NDES following the capture of H~ on the deuteron, over the FSI region, for a 
series of interactions in which the strength of the coupling between the A-A and 3-N is modified 
by multiplying by a factor R, i.e. Ca= — > -RCah with R = 0.5,0.75, 1.0, and 1.25. 
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